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Abstract 

The two component model of coexisting local electron pairs and itinerant fermions 
coupled via charge exchange mechanism, which mutually induces superconductivity 
in both subsystems, is discussed. The cases of isotropic s-wave and anisotropic pair- 
ing of extended s and d x 2_ y 2 -wave symmetries are analyzed for a 2D square lattice 
within the BCS-mean field approximation and the Kosterlitz-Thouless theory. We 
determined the phase diagrams and superconducting characteristics as a function 
of the position of the local pair (LP) level and the total electron concentration. The 
model exhibits several types of interesting crossovers from BCS like behavior to 
that of LP's. Some of our results are discussed in connection with a two-component 
scenario of preformed pairs and unpaired electrons for exotic superconductors. 
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A system of interacting charged bosons (bound electron pairs) and electrons 
can show features which are intermediate between those of local pair super- 
conductors and those of classical BCS systems. Such a two component (boson- 
fermion) model is of relevance for high temperature superconductors (HTS) 
and other exotic superconductors [1,2,3,4,5,6,7,8,9,10]. Recently, we have stud- 
ied a generalization of this model to the case of anisotropic pairing [9,10]. 
Here, we shortly outline the study and present further results concerning the 
phase diagrams and superconducting properties of such a system in the case of 
isotropic s-wave pairing as well as for anisotropic pairings of d- and extended 
s-wave symmetries, for a 2D square lattice. 
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1 The Model 



We consider the model of coexisting electron pairs (formed by "d" electrons) 
and itinerant "c" electrons defined by the following effective Hamiltonian 

h = E( e k - n)<L<** + 2 E( A o - - ]r Jtjbtbj + 

kcr i ij 

E [^q( k ) C k+ q /2,T C -k+ q /2,| 6 q + H - C ] > C 1 ) 
k,q 

where e k refers to the band energy of the c-electrons, A measures the relative 
position of the LP level with respect to the bottom of the c-electron band, \i is 
the chemical potential which ensures that the total number of particles in the 
system is constant, i.e. n = ^(Dko- < c L- c ko- > +2J2i < bjh >) = n c + 2n B - 
n c is the concentration of c-electrons, % is the number of d-pairs (hard-core 
bosons) per site. J^- is the pair hopping integral The operators for local pairs 
b\, bi obey the Pauli spin 1/2 commutation rules. V q (k) describes the coupling 
between the two subsystems. Furthermore, we set = 0. We will analyze the 
case Vq(k) = Vo(k) = I<p k /y/~N, and neglect its q dependence at small q. The 
interaction term takes the form of coupling, via the center of mass momenta 
q, of the singlet pair of c-electrons 5 q and the hard-core boson 6 q : 



#i = ^E 7 (£ q & q + & q £ q )- (2) 

ViV q 

Bq = 0fc c k+ q /2,t c -k+ q /2,| denotes the singlet pair creation operator of c- 
electrons and I is the coupling constant. The pairing symmetry, on a 2D 
square lattice, is determined by the form of <pk, which is 1 for on-site pairing 
(s), <pk = Ik = cos(k x a) + cos(k y a) for extended s-wave (s*) and <pk = Vk = 
cos(k x a) — cos(kya) for d^-^-wave pairing (d). In general, one can consider a 
decomposition I fa = g + g s ^k + gdVk, with appropriate coupling parameters 
for different symmetry channels. 



The numerical results presented below are based on the BCS-Mean-Field Ap- 
proximation (MFA) and the Kosterlitz-Thouless (KT) theory for 2D super- 
fluid. In the BCS-MFA the free energy of the system is calculated to be: 



F/N = ~JN ^ ln [2 cosh (^/2)] - ^ In [2 cosh(/3A)] + C, (3) 

C = -e b + A + fi(n c + 2n B ) - 2/2 - 2I\x \p x , (4) 

where the quasiparticle energy of the c-electron subsystem is given by 

£ k = y/% + Al and e k = e k -/i, A 2 , = / 2 0i(pg) 2 , A = ^(A - ^ + P\x \* , 
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(5 = 1/ksT. The c-electron dispersion is ek = ek— e& = — 2t [cos(/c x a) + cos^a)] — 
4t 2 cos^a) cos(kyd) — 65, with the nn and nnn hopping parameters t and t 2 , 
respectively, e& = mine^ . It should be noted that the energy gap in the c-band 
is due to nonzero Bose condensate amplitude (| < b > \ ^ 0), and well defined 
Bogolyubov quasiparticles exist in the superconducting phase. The energy 
spectrum of the system is characterized by and A. The superconducting 
order parameters: x = jf J2k 4>k < <4t c -A:| > an d Po = < bj + h > and 

the chemical potential \i are given by 



OF OF OF 

The superfluid density derived within the linear response method and BCS 
theory is of the form: 




depth. 

Let us also consider the pair propagator for c-electron subsystem 



G 2 (q,u;) = -i«B q |flt » w> (7) 

where << -B q |-B q >> w is the time Fourier transform of the Green's function: 
-iQ(t - *') < [B q (t),Bj(t')] >• In the normal state, using the equation of 
motion technique and Random Phase Approximation one gets: 



G 2 (q,w) = - T -— ( r, 8 



where 



k a; - efc+ g /2 - 

I 2 

Veff{uj) = cu-2(A -/i) tanh [(3{A ° ~ ^ )] ' (10) 

x(q, a;) is the pair susceptibility and v e ff (u)-the effective interaction mediated 
by LP's. Introducing the generalized T-matrix r(q, u) via the relation G2 = 

X^/veff, one has 
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r(q,c) = - Y u \ . (11) 

An instability of the normal phase occcurs at r _1 (0, 0) = (the Thouless 
criterion), and T c is given by 



1 = 



2(A -//) 



tanh 



c 



ek/2) 



2e k 



(12) 



in full agreement with the MFA expression which one obtains from Eqs.(5), if 
xo -> 0,pg -> 0. 



The MFA transition temperature (Tj^^) at which the gap amplitude vanishes 
yields an estimation of the c-electron pair formation temperature [1]. Due to 
the fluctuation effects the superconducting phase transition will occur at a 
critical temperature lower than that given by the BCS-MFA theory. In 2D, 
T c can be derived within the KT theory for superfluid [11], which describes 
the transition in terms of vortex-antivortex pair unbinding. We evaluate T c 
using the KT relation for the universal jump of the superfluid density p s at T c 
[11,12]: 

-k B T c = p a (T c ), (13) 

71 



where p s {T) is given by Eq.(6) and xq(T), Pq(T), h(T) are given by Eqs.(5). 
Hence, the critical temperature denoted further by T C XT is determined from 
the set of four self-consistent equations. 



2 Results and discussion 



We have performed a comprehensive analysis of the phase diagrams and super- 
fluid properties of the model (1) for different pairing symmetries [10]. Below 
we will discuss these results including the supplementary ones. 
In Fig.l we show the ground state diagrams as a function of n and Aq/D, 
plotted for |J | = (I — — \I \) and two fixed values of t 2 /t. A weak intersub- 
system coupling |/ | will not change much characteristic lines of the diagrams, 
while increasing |7 | expands the range of the LP+E regime [13]. Depending 
on the relative concentration of "c" and "d" electrons we distinguish three 
essentially different physical situations (in the absence of interactions). For 
n < 2 it will be: 

(i) A < such that at T = OK all the available electrons form local pairs (of 
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"d" electrons) (the "local pair" regime, 2n B ^> n c ) (LP); 

(ii) A > such that the "c" electron band is filled up to the Fermi level 
H = A and the remaining electrons are in the form of local pairs (the c+d 
regime or Mixed, < 2n B ,n c < 2) (LP+E); 

(iii) A > such that the Fermi level /i < A and consequently at T = OK all 
the available electrons occupy the "c" electron states (the c- regime or"BCS", 
n c > 2n B ) (E). 

For |7 | 7^ 0, in the case (ii) superconductivity is due to the interchange be- 
tween local pairs of "d" electrons and pairs of "c" electrons. In this process "c" 
electrons become "polarized" into Cooper pairs and "d" electron pairs increase 
their mobility by decaying into "c" electron pairs. To this intermediate case 
neither the standard BCS picture nor the picture of local pairs applies and 
superconductivity has a "mixed" character. The system shows features which 
are intermediate between the BCS and preformed local pair regime. This re- 
gards the energy gap in the single-electron excitation spectrum (E g (0)), the 
k B T c / E g (0) ratio, the critical fields, the Ginzburg ratio /t, the width of the crit- 
ical regime as well as the normal state properties. In case (i) the local pairs 
of "d" electrons can move via a mechanism of virtual excitations into empty 
c-electrons states. Such a mechanism gives rise to the long range hopping of 
pairs of "d" electrons (in analogy to the RKKY interaction for s-d mechanism 
in the magnetic equivalent). The superconducting properties are analogous to 
those of a pure local pair superconductor [1,6]. In case (iii), on the contrary, 
we find a situation which is similar to the BCS case: pairs of "c" electrons 
with opposite momenta and spins are exchanged via virtual transitions into 
local pair states. 

The generic phase diagrams for different pairing symmetries plotted as a func- 
tion of the position of the LP level A at fixed n are shown in Fig.2a-c. In 
all the cases one observes a drop in the superfluid stiffness (and in the KT 
transition temperature) when the bosonic level reaches the bottom of the c- 
electron band and the system approaches the LP limit. In the opposite, BCS 
like limit, T^ T approaches asymptotically T c MFj4 , with a narrow fluctuation 
regime. Between the KT and MFA temperatures, the phase fluctuation effects 
are important. In this regime a pseudogap in the c-electron spectrum will de- 
velop and the normal state of LP and itinerant fermions can exhibit non- Fermi 
liquid properties [5]. 

The analysis of the Mixed-LP crossover indicates that when the LP level is 
lowered and reaches the bottom of the fermionic band an effective attraction 
between fermions becomes strong, since it varies as J 2 /(2A — 2ji) and fi ~ A 
(Eq. 10). In this regime the density of c electrons is low and formation of bound 
c-electron pairs occurs. It gives rise to an energy gap in the single-electron 
spectrum independently of the pairing symmetry. We evaluated the binding 
energies of c- electron pairs and found that T c Mi?A scales with the half of their 
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binding energy for Ao < 0. The superconducting transition temperature is 
here always much lower than the c-pair formation temperature (T C MFA ) and 
rapidly decreases with |A /D|. In such the superconducting state can 

be formed by two types of coexisting (composite) bosons: preformed c- electron 
pairs and LP's [10,14]. 

Comparing T c vs Ao plots for various pairings one finds that in the case of nn 
hopping only, the d and s -wave pairings are favorable for higher concentration 
of c-electrons, while the s*-wave can be stable at low n c . This is clearly seen 
from Fig.3, where the plots of T c Mi?A vs A for s, s* and rf-pairings are given 
for n — 1, \Iq\/D = 0.25, and from Fig. 4 which shows variation of n c and % 
with Aq/D for the same values of n and Aq/D for d -wave pairing (for s and 
s* -wave symmetries the plots of n c and % are very similar). 

A strongly nonmonotonous behavior of T c vs A for s*-wave pairing and n > 1 
(cf. Fig. 2b) reflects the fact that this type of pairing in a single band model 
is realized for low electron (or hole densities), if ti = [1]. In general, with 
increasing n c (A /D) one can have a sequence of transitions: s* — > d (or 
s)— > s* (if 1 < n < 2), or the transitions: s* — > d (or s) (if n < 1). 
The nnn hopping t 2 (with opposite sign to t) can strongly enhance T c for d- 
wave symmetry, moreover it favors the d and s-wave pairings for lower values 
of n c [10]. 

The region between T c Mi?A and T^ T , where the system can exhibit a pseu- 
dogap, expands with increasing intersubsystem coupling |J |. It is seen from 
Fig. 5 that, except for \Iq\/D <C 1 the coupling dependences of T^ FA and 
T^ T are qualitatively different. T c Mi?A is an increasing function of |io| for all 
the pairing symmetries. On the other hand, T^ T vs \Iq\ increases first, goes 
through a round maximum and then decreases (similarly as it is observed in 
the attractive Hubbard model [12]). The position of the maximum corresponds 
to the intermediate values of \Iq\/D and it depends on the pairing symmetry 
as well as the values of Aq/D and n. For large |/o|, the T^ T are close to the 
upper bound for the phase ordering temperature which is given by 7rp s (0)/2. 

As for the evolution of superconducting properties with increasing n one finds 
three possible types of density driven changeovers [10,13]: (i) for 2 > A /D > 
0, "BCS" — ► Mixed — ► "BCS"; (ii) for A /D > 2: "BCS" — ► "LP" and (iii) 
for A /D < 0: "LP" — > "BCS". Only if the LP level is deeply located below 
the bottom of c-band, the system remains in the LP regime for any n < 2. 

Some of our findings can be qualitatively related to experimental results for 
the cuprate HTS where a pseudogap exists. It has been suggested by ARPES 
experiments, that for underdoped cuprates the Fermi surface in the pseudo- 
gap phase is truncated around the corners due to the formation of preformed 
(bosonic) pairs with charge 2e, whereas the "electrons" on the diagonals re- 
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main unpaired [7a]. In the present two component model such a situation is 
obtained when LP's and c-electrons coexist in the mixed regime. The linear 
T-dependence of the superfluid density has been observed experimentally in 
copper oxides and also in several organic superconductors. This points to an or- 
der parameter of ci^-j^-wave symmetry and existence of nodal quasiparticles. 
In the present model the gap ratio is nonuniversal and can deviate strongly 
from BCS predictions (particularly in the d-wave case for which it is always 
enhanced) [10]. This feature is also found for several exotic superconductors. 
The Uemura plots i.e., the T c vs zero-temperature phase stiffness p s (0) (and 
the scaling T c ~ p s (0)) reported for cuprates and organic superconductors can 
be reproduced within the model for extended s- and <i-wave order parameter 
symmetry ([9], Figs. 11, 13 in Ref. [10]). 

Finally, we point out that to relate more quantitatively the present two- 
component scenario to the phase diagram of cuprates one has to consider 
doping dependences of the model parameters. 
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Fig. 1. Ground state diagram of the model (1) as a function of n = n c + 2n# versus 
Aq/D (D = zt = 4t) for | Jo | = 0, plotted for a square lattice with t% = (thick 
solid lines), rectangular DOS (dashed lines) and with t%jt = —0.45 (thin lines). LP - 
nonmetallic state of LP's (2n# = n, n c = for n < 2, 2ns = 2— n, n c = 2 for n > 2), 
E - metallic state of c -electrons (n c = n, 2n£ = for n < 2, n c = 2 — n, 2ub = 2 
for n > 2), LP+E - mixed regime (0 < n c , 2ub < 2). 
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Fig. 2. Phase diagrams of the induced pairing model as a function of Aq/D 
derived for (a) s-wave, (b)- extended s-wave, (c) ci-wave symmetry (£2 = 0). 
In all the figures the dashed line shows the BCS-MFA transition tempera- 
ture, while the line with diamonds the KT transition temperature calculated for 
n = 1.5,\Iq\/D = 0.25. LPN-normal state of predominantly LP's, EM-electronic 
metal, LPS+ES-superconducting (SC) state, PG - pseudogap region. 
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Fig. 3. T^ F as a function Aq/D for different pairing symmetries plotted for 
n = 1,\Iq\/D = 0.25, ti = 0. Solid lines are the results for the d = 2 square 
lattice, while the dashed ones - for the quasi 2D case with tj_/t = 0.1 
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Fig. 4. (a) Variation of n c (concentration of c-electrons) , ns (concentration of LP's) 
and superconducting order parameters Pq,xo at T = as a function of Aq/D, for 
n = 1, |/ol/-^ = 0.25, t2 = 0. d-wave pairing, (b) The corresponding plots of n c ,n^ 
and n/D at T^ FA . 
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Fig. 5. Temperatures J 1 * and T c as a function of \Iq\/D for various types of 
pairing for a 2D lattice (£2 = 0). The lines with full symbols and thick solid line 
denote the KT transition temperatures. The line with empty squares is T^ FA for 
d-wave pairing (Aq/D = l,n = 2), while the line with open circles is T^ FA for 
s*-wave pairing (Aq/D = 0.5, n = 1). 
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